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SYMBOLS 


C L 

D 

d 

d c 


dcjd 

dj 


lift coefficient 
drag force 

horizontal distance subtended by an altitude 
segment 

horizontal distance subtended by the constant 
speed leg of a segment in the speed profile (see 
fig. 8) 

horizontal distance subtended while climbing or 
descending in an altitude segment 

horizontal distance subtended by the final leg 
of a segment in the speed profile (see fig. 8) 


d 0 horizontal distance subtended by the first leg 

of a segment in the speed profile (see fig. 8) 


di 

E 

EPRmax 
EP Rmin 

k 

K 

f 

h 


h 


the horizontal distance subtended while flying 
along a segment of the speed profile 

horizontal distance subtended during level flight 
in an altitude segment 

total energy of aircraft 
maximum EPR setting 
minimum EPR setting 
aircraft energy rate 
normalized energy rate 
fuel consumption 

fuel consumption for the first leg of speed seg- 
ment 

fuel consumption for the last leg (usually the 
third leg) of speed segment 


fuel consumption for the constant-indicated- 
airspeed leg of speed segment 

fuel consumption for the constant-true-airspeed 
leg of speed segment 

fuel consumption rate of a speed segment 
acceleration of gravity 

altitude at the end of the first leg of a speed 
segment 

altitude at the beginning of the last leg of a 
speed segment 

altitude at the end of a speed segment 

altitude at the the beginning of the i th altitude 
waypoint 

altitude at the the end of the i th altitude 
waypoint 

altitude at the beginning of a speed seg- 
ment 

current aircraft heading 

final heading for standard route or runway head- 
ing 

reference aircraft heading at waypoint i 
index for altitude segment pattern (see table 2) 

mass of the aircraft 
Mach number 

number of horizontal waypoints in a horizontal 
flightpath 

number of altitude waypoints in an altitude 
profile 

number of speed waypoints in a speed profile 

transition point from straight flight to turn for 
the horizontal flightpath segment i 


jO S3 ft) 


P 


Qi 


S 

(Sj , hj ) 

(«*,VJk) 

S 

T 

Tm ax 

Tmin 

tci 

let 

td 

tmax 

tmin 

tnom 

toa 

h 


transition point from turn to straight flight for 
the horizontal-flightpath segment i 

turning radius 

turning radius for waypoint i 
turning radius for the last waypoint 

distance along the horizontal flightpath, a vari- 
able 

coordinates of altitude waypoint, the first com- 
ponent is the distance from touchdown and the 
second component is the altitude 

coordinates of speed waypoint, the first com- 
ponent is the distance from the touchdown 
point, the second component is the indicated 
airspeed 

wing reference area, 1,560 ft 2 for the example 

aircraft 

thrust 

maximum thrust 
minimum thrust 

time duration for traversing the constant- 
indicated-airspeed leg of a speed segment 

time duration for traversing the constant-true- 
airspeed leg of a speed segment 

the desired time of arrival 

maximum time duration for flying along a speed 
segment 

minimum time duration for flying along a speed 
segment 

time duration for flying along a speed seg- 
ment 

the controlled time of arrival 

time duration for traversing the first (nondegenerative) 
leg of a speed segment 


IX 


t& time duration for traversing the last (usually the 

third) leg of a speed segment 

V a aircraft velocity, a variable 

V c i true airspeed at the beginning of a constant- 

indicated-airspeed leg of a speed segment; its 
corresponding indicated airspeed is denoted by 

v. 

V C 2 true airspeed at the end of a constant-indicated- 

airspeed leg of a speed segment; its cor- 
responding indicated airspeed is denoted by 
V'c 

Vf final true airspeed in a segment of a speed 

profile 

V g ground speed of aircraft 

Vi as indicated airspeed 

V 0 initial true airspeed in a segment of a speed 

profile 

V oa speed for the constant leg in a controlled-time- 

of-arrival speed segment 

V w wind velocity 

V a dVa/dt 

V w dV w / dt 

V' c constant indicated airspeed in a segment of 

speed profile (see fig. 8) 

Vj final indicated airspeed in a segment of a speed 

profile 

V k the initial indicated airspeed in the k th speed 

segment 

V 0 initial indicated airspeed in a segment of a speed 

profile 

W aircraft weight 

Xa>ya x,y coordinate of aircraft position, runway 

reference 


X 


Xi,Vi 

xpi,ypi 
X(Ji, VQi 
S f 

AC D,Ssp 

Afi 


A Si 
A ti 

e t,k 

e 

7 

7 

Iclimb 

^descent 

Idz 


7 ez 


p(h) 

Pel 

a 

tymax 


x, y coordinates of waypoint i, runway reference 
x,y coordinates of point P*, runway reference 
x, y coordinates of point Q t , runway reference 
flap deflection, deg 

incremental change in drag coefficient due to 
employment of speed brake 

incremental fuel consumption for i th step of in- 
tegration in a deceleration/ acceleraton leg of a 
speed segment 

incremental distance or step size for integrating 

speed leg 

incremental time for integrating speed leg 

the allowable time error for the k th segment of 
a controlled-time-of- arrival profile 

atmospheric temperature; used for computing 
TSFC 

ratio of specific heat = 1.4 
flightpath angle 

specified flightpath angle for climb 
specified flightpath angle for descent 

angle of depression, arc tangent of altitude 
drop per unit horizontal distance in the altitude 
profile 

angle of elevation, arc tangent of altitude rise 
per unit horizonal distance in the altitude 
profile 

atmospheric density as a function of altitude 
atmospheric density at sea level 
speed of sound 
maximum bank angle 


xi 


Subscripts: 
a aircraft 

i waypoint i for horizontal flightpath 
t i th integration step 
j waypoint j for altitude profile 
k waypoint k for speed profile 
n the last horizontal waypoint in a route 

Abbreviations and Acronyms: 

ATC air traffic control 

CAS calibrated airspeed 

C/D climb/descent in the altitude profile 

CLA cleared for approach 

EPR engine pressure ratio 

IAS indicated airspeed 

KIAS indicated airspeed in knots 

ILS instrument landing system 

LF level flight in the altitude profile 

OM outer marker 

ST straight flight followed by a turn 

STS straight flight, turn and then straight flight 

TAS true airspeed 

TD touchdown or touchdown point 

TOA time of arrival 

TS turn followed by straight flight 

TSFC thrust specific fuel consumption 

TST turn-straight flight and then turn 

WP waypoint 

4D four-dimensional guidance, x,y,z coordinates 

and time 
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SUMMARY 


This report describes concepts and computer algorithms for generat- 
ing time-controlled four-dimensional descent trajectories. The algorithms 
have been implemented in the NASA Ames Air Traffic Control Simulator 
and used by experienced controllers in studies of advanced air-traffic- 
flow management procedures. A time-controlled descent trajectory com- 
prises a vector function of time, including position, altitude, airspeed, 
and heading, that starts at the initial position of the aircraft and ends 
at touchdown. The trajectory provides a four-dimensional (4D) reference 
path which will cause an aircraft tracking it to touchdown at a predeter- 
mined time with a minimum of fuel consumption. The problem of 
constructing such trajectories is divided into three subproblems involv- 
ing synthesis of horizontal, vertical, and speed profiles. The horizon- 
tal profile is constructed as a sequence of turns and straight lines pass- 
ing through a specified set of waypoints. The vertical profile con- 
sists of a sequence of level-flight and constant-descent-angle segments 
defined by altitude waypoints. The speed profile is synthesized as 
a sequence of constant-Mach-number, constant-indicated-airspeed, and 
acceleration/deceleration legs. It is generated by integrating point-mass 
differential equations of motion, which include the thrust and drag models 
of the aircraft. A speed profile with a specified touchdown time along a 
known horizontal profile is obtained by iteration on the constant-Mach- 
number and indicated-airspeed segments of the descent. Three types of 
4D trajectories — referred to as standard approach, 4D capture, and 
cleared-for-approach — are generated by different versions of the al- 
gorithm. 

The algorithm also distinguishes between trajectories for 4D-flight- 
management-equipped and unequipped aircraft. Trajectories for equipped 
aircraft have more nearly fuel-optimal altitude profiles and are assumed to 
fly uninterrupted from starting point to touchdown. Trajectories for un- 
equipped aircraft have altitude profiles consistent with current operational 
practice. 


INTRODUCTION 


The air-traffic-control (ATC) system is in the midst of changes in 
systems and procedures to cope with increasing air-traffic congestion at 
hub airports and to achieve energy conservation. An on-board system 
has been developed by Ames Research Center which can assist with both 
fuel and capacity problems. The system, which is specified by a set of 
algorithms to be programmed for use on board the aircraft, generates 
a four-dimensional (4D) trajectory that optimizes fuel use along an as- 
signed route and meets an assigned time requirement. The ability to 
arrive at an assigned time permits a more orderly, fuel-efficient flow 
of traffic than does the standard vectoring mode. Furthermore, if 
a sufficient number of 4D-equipped aircraft are present in a terminal 
area, the orderly and predictable traffic flow can be used to improve 
capacity. 

The 4D aspect of the on-board system evolved from conceptual 
studies through simulation and flight-test experiments. The 4D guidance 
and control concept was described in references 1 and 2, and was later 
simplified to make it suitable for real-time simulation (ref. 3). The tech- 
nique was implemented in an airborne computer and flight tested on a low- 
performance aircraft in the terminal area (ref. 4). It has also been imple- 
mented for air-traffic-control simulation studies using experienced ATC con- 
trollers (ref. 5). In this report, previously developed techniques (refs. 1-8) 
are combined with new 4D algorithms that span the full operational en- 
velope of a commercial jet transport, from the top of descent to touch- 
down. 


Briefly, the on-board system operates as follows: The individual 
aircraft enters the terminal area at the feeder-fix and captures one of 
the standard terminal-area routes. Capturing means a smooth transi- 
tion from current aircraft state (which includes the x, y, z position coor- 
dinates, heading, and speed) to some specified state along the route. All 
routes are 4D-guidance routes with a time slot for landing assigned to 
individual aircraft. Next, the aircraft flies along the route until touch- 
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down, unless the controller interferes by vectoring the aircraft. The 
vectored aircraft can be reassigned to a 4D route via the capture com- 
mand. 


This report provides a detailed description of the 4D algorithm as 
implemented in the real-time ATC simulation at Ames Research Center. 
The report is organized as follows: The first section defines a 4D trajectory 
and describes the algorithm for synthesizing a standard 4D, terminal- area 
route. The algorithm contains two parts, one part for synthesizing the 
three-dimensional (3D) path composed of horizontal and altitude profiles 
and another part for synthesizing speed profiles for controlled time of 
arrival. The next section describes the modification to the algorithm 
for synthesizing trajectories for the cleared-for-approach (CLA) command. 
This command causes an aircraft to capture the instrument-landing-system 
(ILS) glide slope and localizer and fly the final approach to touchdown 
(TD). 


SYNTHESIZING A 4D TRAJECTORY 


The components of a 4D trajectory are the horizontal flightpath, 
the altitude profile, and the speed profile. For the purpose of synthesis, 
they are essentially independent of each other, but there is some coupling 
between them. The horizontal distance subtended by the horizontal 
flightpath is used as the independent variable for specifying the altitude 
and speed profile. The altitude profile and speed profile are coupled 
by the energy-rate capability of the aircraft. The turning radius used 
for specifying the horizontal flightpath depends on the true airspeed 
(TAS) and the wind speed. The speed profile is usually specified in 
indicated airspeed (IAS), and the conversion from IAS to TAS is altitude- 
dependent. 

The horizontal flightpath is specified by the initial route heading, 
the final route heading, the ( x , y) coordinates of the waypoints, and the 
turning radius at each waypoint. The altitude profile is specified by 
selecting altitude waypoints along the horizonal flightpath, and choos- 
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ing the altitude and the minimum and maximum flightpath angle at 
these points. Similarly, the speed profile is specified by selecting speed 
waypoints along the horizontal flightpath, and choosing the maximum 
and minimum speed at these points, as well as the desired arrival time 
at touchdown. From these specifications, the synthesis procedure fol- 
lows. First, synthesize the horizontal flightpath and then locate the al- 
titude and speed waypoints along the horizontal flightpath. Second, syn- 
thesize the altitude profile; and third, use the altitude profile to deter- 
mine the forces, namely thrust, lift, and drag, to synthesize the speed 
profile for the specified time of arrival. After the speed profile is syn- 
thesized, the fuel consumption for the profile is known. Since the turn- 
ing radii for specifying the horizontal flightpath are speed-dependent, their 
values will be checked to make certain that the turns specified are flyable, 
which is usually the case. If they are not flyable, some iteration is neces- 
sary. 

The following definitions are adopted in the development of 
the trajectory synthesis methodology. A profile denotes the entire 
airspeed or altitude variation as a function of distance starting at the 
aircraft position and ending at touchdown. A segment is the por- 
tion of the profile or flightpath between two adjacent waypoints. A 
leg is a further subdivision of a segment. A detailed discussion of 
the synthesis technique for the horizontal flightpath and the altitude 
profile follows. The speed profile will be discussed in a later sec- 
tion. 

3D Flightpath Patterns 

The horizonal flightpath is formed by a concatenation of horizontal 
flightpath segments. There are three types of segments, straight-flight/turn 
(ST), straight- flight/turn/ straight-flight (STS), and turn/straight- flight/turn 
(TST). The altitude profile is formed from a concatenation of 
up to four altitude segments laid out along the previously defined 
horizontal segments. 

Figure 1(a) shows the pattern for the ST segment. It is defined by a 
pair of waypoint coordinates (x n —i,y n — 1 ) and ( x n ,y n ), the final heading 
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Hj , and a turning radius R. The turning radius is a function of V g , the 
ground speed of the aircraft, and <t> max the maximum bank angle of the 
aircraft; that is, 


R = 


V 2 

v 9 

g tan <f>max 


( 1 ) 


where g is the acceleration of gravity. The segment begins with waypoint 
n-1 and is completed at waypoint n, with the aircraft heading lined up 
with Hf. The point P n is the transition point from straight flight to 
turn. Its coordinates are denoted by (xp n —i,yp n — i). Figure 1(b) shows 
the pattern for the STS segment. It is specified by the coordinates of 
three waypoints, (21,3/1), (x 2 ,y 2 ), and (^3,2/3), and the turning radius 
R, which is computed by using equation (1). The pattern starts at 
waypoint 1, approaches but does not pass through waypoint 2, and finishes 
at waypoint 3 . The point Q is another transition point from turn to straight 
flight. Figure 1(c) shows the pattern for the TST segment. It is specified 
by the coordinates of two waypoints, (21,3/1) and (22,2/2), the headings at 
these two waypoints, Hi and H 2 , and the turning radii at these points, R\ 
and R 2 • 


Synthesizing Horizontal Flightpath Segments 

Standard route- The complete horizontal flightpath for a stan- 
dard route is specified by the initial heading of a selected route Hi, 
the corresponding initial position (21,3/1) on that route, referred to as 
the feeder- fix; the runway heading Hf) the coordinates of the touch- 
down point (2 n , 2 /n); and the coordinates of intermediate waypoints, 
(22, 2/2)— (^n— 1 , Vn 1 )• It is synthesized backward, starting from the touch- 

down point, through intermediate waypoints and finishing at the initial 
position. The first segment is synthesized as an ST pattern in the for- 
ward direction or as turn/straight-flight (TS) pattern in the backward direc- 
tion from known (2 n ,y n ), (z n _ \,y n — 1), Hf and R n , the turning radius 
specified for waypoint n. After the segment is synthesized, the heading H n , 
formed by the line between {x n —i,y n —i) and the point P n (see fig. 1(a)) 
is known, where P n is the transition point from ST (or the point for the 
beginning of a turn) in the forward direction. Its coordinates, denoted by 
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Figure 1.— Patterns for horizontal flightpath segments. 


(xp n ,yPn), are also known. The intermediate segments are synthesized 
backward either as an STS or a TS segment. The first segment (the 
last segment in the backward direction) is always synthesized as a TST 
pattern, using (21,3/1), (22,3/2), #1 , and # 2 as input. This segment is 
essentially treated as a capture route and is explained in the next sec- 
tion. 
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Selection of an STS or a TS pattern for synthesizing intermediate 
segments is based on the feasibility of a solution and the amount of 
heading change required at a waypoint. In general, if the included 
angle, the angle formed by the two line segments connecting a waypoint 
to its two adjacent neighbors, is less than 90°, a TS pattern is al- 
ways used. If the included angle is greater than 90°, in which case 
the heading change is less than 90°, the procedure is to attempt syn- 
thesis of an STS pattern. If there is no solution, a TS pattern is 
also used. A solution may not exist for the STS pattern because turns 
at adjacent waypoints may overlap. This problem is explained further 
in a later section (Cleared-for-Approach Command). References 3 and 
6 derive the synthesis relation for the STS and TS patterns, respec- 
tively. 

Capture Route- This algorithm solves the problem of synthesiz- 
ing a 4D trajectory for an aircraft not starting at a standard feeder- 
fix location. The uncertainties of weather, avoidance of other traffic, 
and holding, among other factors, can result in an aircraft deviat- 
ing from its standard route, both before and after departing the 
feeder-fix. The capture algorithm generates a reference trajectory 
to return the aircraft to a standard route simply by specifying a 
capture waypoint on that route. It is also used for generating 
the segment between the first and second waypoints of a standard 
route. 


Since the algorithm is basically identical to that described previously 
in reference 6, it will only be summarized here. The technique 
developed in reference 6 addresses the problem of transferring an aircraft 
from a given initial position and heading to some final position and 
heading. There exist at least two and at most four solutions to 

this problem. The general solution yields a TST pattern. The al- 
gorithm resulting from application of this technique finds all the solu- 
tions first and then selects the one that yields the minimum path 
length. 


After the horizontal profile has been synthesized to the capture 
waypoint, the non-standard initial position may be considered as defining 


7 


the start of another route. The ST pattern is also synthesized by 
this algorithm as a degenerative case, by setting the appropriate turn- 
ing radius to zero and consequently ignoring the corresponding head- 
ing. 


Examples of synthesized routes - Shown in figure 1(d) are the 
different horizontal paths produced by the two synthesis methods, STS and 
TS. The STS method rounds the corner at a waypoint, so the waypoint 
itself is not on the path. On the other hand, the TS pattern creates a path 
that passes through the waypoint. 


wp 3 



Figure 1.— Concluded. 
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Figure 2 shows a synthesized route for an airport similar to 
that at the Kennedy Airport in New York. Table 1 shows the cor- 
responding numerical values for the route. Since air-traffic controllers 
and computer programmers are interested in different aspects of the al- 
gorithm for ATC simulation, two sets of numbers are used to desig- 
nate the waypoints. The first set refers to the numbers as shown in 
the terminal-area map (for the benefit of the air-traffic controller); it 
is shown in column 1. The second set is arranged in numerical order 
and refers to elements stored in an array in computer programs; it is 
shown in column 2. For example, waypoints 54 and 39 are, respec- 
tively, the first and last waypoints of Ellis route as shown in column 1. 



Figure 2.— Synthesized horizontal flightpath. 
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TABLE 1 -SYNTHESIZED HORIZONTAL FLIGHTPATH 


Way- 

point 

No. 

Waypoint 
coordinates, 
n. mi. 

Heading, 

deg 

Straight- 

flight 

distance. 

Turn 
distance, 
n. mi. 

Distance 
to TD, 
n. mi. 

Remarks 






n. mi. 




X 

y 





54 

1 

84.45 

163.70 

120.4 

9.27 

0 

118.96 

ST4 Ellis approach 

45 

2 

92.44 

159 

120.4 

15.57 

.56 

109.69 


35 

3 

106.13 

151 

137.9 

27.60 

.28 

93.56 

SAX Sparta VOR 

22 

4 

124.92 

130.22 

128.9 

26.75 

2.92 

65.68 

Ellis intersection 

19 

5 

147.31 

112.13 

-138.5 

11.45 

2.77 

36.01 


17 

6 

139.92 

100.84 

-49.8 

7.47 

2.04 

21.79 


9 

7 

129.32 

107.43 

14 

3.79 

.86 

12.28 


6 

8 

130.62 

112.63 

41.4 

4.53 

.05 

7.63 

1 LS Intercept 

3 

9 

133.91 

116.38 

43 

3.05 

0 

3.05 

Outer marker 

39 

10 

136.01 

118.62 

43 



0 

TD point at runway 4R 


The value of the x,y coordinates will be designated by {xi,yi) and 
(xio,Vio)l ^e subscripts are shown in column 2. Columns 3 and 4 
show, respectively, the x, y coordinates of the waypoints in nautical miles; 
column 5 the heading at the waypoint in degrees, columns 6 and 7 
the straightflight distance and the distance subtended while making the 
specified turn respectively; column 8 shows the distance to go from touch- 
down. 


The route is specified by the (x, y) coordinates of n waypoints, the 
initial heading Hi , and the final (or runway) heading H / . The numerical 
value of n is 10, and the numerical values of Hi and Hj are 120° and 43°, 
respectively. 

The last segment of the route is synthesized forward in the backward 
direction (i.e., last segment first) as a TS pattern. It is specified by Hj, 
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and the coordinates of waypoints 39 and 3, that is ( 210 , 2 / 10 ) and ( 29 , 2 / 0 ), 
where 

( 210 , 2/10 ) = (136.01, 118.62) (2a) 

( 29 , 1 / 9 ) = (133.91, 116.38) (26) 

The resulting pattern is to fly straight 3.05 n. mi. and turn 0 n. mi., a 
degenerative case of TS pattern consisting of a straight-flight leg only. Since 
the turn is zero, the transition point P 10 coincides with waypoint 39, the 
touchdown point. The next segment is synthesized as a STS segment since 
the included angle is greater than 90° (in fact almost 180°). It is specified 
by the coordinates of waypoints 8 and 9 and of transition-point 10, that 
is, 


Pio,mo) = (136.01,118.62) 

(3d) 

( 20 , 2 / 0 ) = (133.91, 116.38) 

(36) 

(2 S , y 8 ) = (130.62, 112.63) 

(3c) 


The synthesis procedure continues backward until the first segment 
is reached. It is synthesized as a TST pattern, the so-called horizontal cap- 
ture flightpath, since both the initial heading, the final heading, and the 
coordinates of the initial and final waypoints for the segment are specified. 
They are 


Hi = 120.4° 

(4») 

H 2 = 120.4° 

(46) 

(* 1 , 2/i ) = (84.45, 163.70) 

(5a) 

(* 2 , 2 / 2 ) = (92.44,159.00) 

(56) 
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Altitude Profile 


An altitude profile is constructed from a concatenation of altitude 
segments. Four typical segments are shown in figures 3(a)-3(d). Each consist 
of two legs, a climb/descent leg followed by a level-flight leg or vice versa. 
A degenerate altitude segment may contain only one leg, a climb, a descent, 
or level flight. 

In the figure 3, h\ and h% are the assigned altitudes at two 
adjacent altitude waypoints, d is the horizontal distance separating 
these points, d c /d is the horizontal distance subtended during al- 
titude change, and di is the horizontal distance subtended during level 
flight. 


An altitude segment is specified by hi , h 2 , d, idimb, or 7 descent , 
and an index //, • There are two possible values for Ih, 0, and 1. 
This index assigns the priority to either making the altitude change 
first and level flight last, or vice versa. The value of /* and its 
corresponding altitude patterns are summarized in table 2 for the climb 
and descent cases. In the descent case, a value of 0 specifies the pattern 
shown in figure 3(a), level flight first and then descent; and a value 
of 1 specifies the pattern shown in figure 3(c), descent first and then 
level flight. Similarly, in the climb case, a value of 0 specifies the 
pattern shown in figure 3(d), level flight and then climb and a value 
of 1 specifies the pattern shown in figure 3(b), climb and then level 
flight. 


Synthesizing Descent Segment- Consider the descent/level-flight 
case as shown in figure 3(c). Define by 7 ^ the angle of depression, which 
is the arc tangent of the specified altitude change per horizontal distance, 
that is 


7 dz = tan 1 — — — (6) 

Comparing this angle with the given angle of descent yields one of three 
results: there is more than enough, just enough, or not enough horizontal 


12 



a) LEVEL FLIGHT-DESCENT SEGMENT 



d 

b) CLIMB-LEVEL FLIGHT SEGMENT 




d c/d 


h 2 


d e 


H 


d 

c) DESCENT-LEVEL FLIGHT SEGMENT 



d) LEVEL-FLIGHT CLIMB SEGMENT 

Figure 3.— Altitude profile. 
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distance for the specified altitude change. These conditions can be formal- 
ized by the relations, 

Idz ^ 7 descent (7a) 

Idz — ^ descent (76) 

' Idz ^ 7 descent (7 c) 

The comparison assumes that both angles are negative. Therefore, 
equation (7a) means that the magnitude of 7 dz is smaller than 7 descent • 
In this case the pattern consists of a descent leg and a level-flight leg. The 
horizontal distance subtended by each of the legs is determined by the rela- 
tion 


h 2 hi 

dc/d — . 

tail *7 descent 

(8) 

di **— d d c jd 

(9) 


If equation (7b) is true, then the pattern degenerates into one descent leg 
and 

d c /d = d (10a) 

di = 0 (106) 

Finally, if equation (7c) is true, the segment also consists of 
one descent leg but there is not sufficient horizontal distance to com- 
plete the altitude change. Therefore, the descent is propagated into 
the next segment. This condition necessitates the readjustment of 
the specified altitude h 2 to reflect the actual altitude to be achieved 
in the descent. Denote by h’ 2 the adjusted altitude (see fig. 4), 
then 

h r 2 = hi + d tan 7 descent (11) 
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The descent will continue to propagate into adjacent segments until equa- 
tion (7a) or equation (7b) is true (since the descent profile is generated from 
lower altitude to higher altitude, it propagates to the previous adjacent seg- 
ments). If not, the synthesis is terminated and a message will appear on the 
controller display stating No Capture — Altitude Profile. Figure 4 shows 
a two-segment descent profile, where the descent continues to the second 
segment. 


Figure 4.— Two segment altitude profile, descent leg propagated to the previous segment. 


Synthesizing Climb Segment- Next, consider the climb segment 
shown in figure 3(b). Similar to the descent case, define by 7 ez the angle 
of elevation, which is the arc tangent of the specified altitude change per 
horizontal distance, that is 


Comparing this angle with the given angle of climb yields one of the three 
results: more than enough, just enough, or not enough horizontal distance 
for the specified altitude change. These conditions can be formalized by the 
relations 



d 



( 6 ’) 
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lez ^ Iclimb 


(12a) 

Hez == Hclimb (126) 

'7 ez ^ r )climb (12c) 


The comparison now assumes that both angles are positive. If equa- 
tion (12a) is true, then the pattern consists of a climb leg and a level- 
flight leg. Horizontal distances subtended by the two legs are, respec- 
tively, 


d c /d = 


h>2. — h\ 

tan 'Iclimb 


( 8 ’) 


and 


— d ^c/<f 


(9’) 


If equation (12b) is true, then the segment degenerates into one climb leg 
only, and 


•c/d — d 

(10a) 

di = 0 

(106) 


Finally, if equation (12c) is true, the segment also consists of one climb 
leg but there is not sufficient horizontal distance to complete the climb. 
Therefore, the climb is propagated into the next segment. This condition 
also necessitates the readjustment of the specified altitude /12 to reflect the 
actual altitude that can be achieved in the climb. Denote by h' 2 the adjusted 
altitude, then 

h' 2 = h x -f d tan 7 c u m b (13) 

The climb will continue to propagate into subsequent segments until 
equation (12a) or (12b) is true. If not, the synthesis is terminated and 
an appropriate message is displayed for the controller. Figure 5 shows 
a two-segment climb profile, where the climb continues into the second 
segment. 
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Figure 5 - Two segment climb profile - climb propagated into the next segment. 


Two degenerative cases are a segment containing only a climb leg 
and only a descent leg. A third case is a segment containing a level-flight 
leg only. It is specified by assigning the same value to the initial and final 
altitude in the segment, that is 


hi = h 2 (14) 

The altitude profile is specified by the integer n a , the number of 
altitude waypoints, ( Sj,hj ), the waypoint coordinates; the angle of 
climb or descent, 7^; and the altitude-change priority index, . 
The subscript j refers to the ordinal waypoint number, and the coor- 
dinate Sj is the distance from touchdown for the waypoint j. A 
positive value of 7 j is for climb and a negative value for des- 
cent. The altitude profile patterns are summarized and explained in 
table 2. 

The four patterns discussed, the three degenerative cases, the 
ability to propagate climb and descent legs to adjacent segments, and 
the freedom to specify altitude waypoints which are completely in- 
dependent of the synthesized horizontal flightpath, provide the flexibility 
for synthesizing altitude profiles. This level of flexibility is useful 
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TABLE 2.- ALTITUDE PATTERN FOR 
ASSIGNED VALUE OF l h 


'h 

Altitude 

change 

Pattern 

Shown in 

0 

Climb 

Level-flight/climb 

Fig. 3(d) 

1 

Climb 

Climb/level-flight 

Fig. 3(b) 

0 

Descent 

Level-flight/descent 

Fig. 3(a) 

1 

Descent 

Descent/level-flight 

Fig. 3(c) 


during the simulation phase in an ATC terminal-area study experi- 
ment, since altitude-profile specification requires only data entries, and no 
software modifications are required to implement changes in the altitude 
profile. 

Examples of a synthesized altitude profile for the previously syn- 
thesized horizontal flightpath are shown in figures 6 and 7. It should 
be emphasized that for the purpose of synthesizing a 4D trajectory, 
the two are essentially independent of each other except that the 
first and the last waypoints of the horizontal flightpath must coin- 
cide with those of the altitude profile. The first waypoint is the 
initial point of the route and the last waypoint is the touchdown 
point. The first waypoint is 118.96 n.mi. from the touchdown point 
(see table 1). Two other horizontal waypoints of interest are the ILS 
intercept point (waypoint 6) at 7.63 n.mi. and the outer marker 
(OM) (waypoint 3) at 3.05 n.mi. from the touchdown point, respec- 
tively. 


Figure 6 shows the synthesized one-segment altitude profile for the 
horizontal flightpath. This type of profile is used for 4D-equipped aircraft 
because it is a good approximation to the fuel-optimum profiles generated 
by on-board flight-management system (ref. 7). 

For unequipped aircraft, profiles are chosen to approximate cur- 
rent operating procedures by adding a level-flight leg between the 
ILS intercept and the outer marker. This ensures that the aircraft 
will capture the glide slope from below after having captured the 
localizer. 
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Figure 6.— Synthesized altitude profile for equipped aircraft. 


Figure 7 shows a synthesized two-segment profile for unequipped 
aircraft. The numerical values for profile specification and for the syn- 
thesized profile are shown in table 3. 



Figure 7.— Synthesized altitude profile for unequipped aircraft. 
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TABLE 3.- ALTITUDE PROFILES FOR UNEQUIPPED AIRCRAFT 


Segment 

no. 

Distance 
to TD, 
n. mi. 

Altitude, 

ft 

Flightpath 

angle, 

deg 

Level-flight 
distance, 
n. mi. 

Climb- 

descent 

distance, 

Pattern 


n. mi. 


Two-segment profile 


118.96 

30,000 





i 

7.63 

993 

-3 

10.81 

100.52 

Level flight, descent 

2 

0 

12 

-3 

4.58 

3.05 

Level flight, descent 

Three-segment profile 


118.96 

33,000 





i 

7.63 

993 

-3 

10.81 

100.52 

Level flight-descent 

2 

3.05 

993 

0 

4.58 

0 

Level flight 

3 

0 

12 

-3 

0 

3.05 

Descent 


The synthesized profile contains four legs, including a level-flight 
segment of length 4.58 n.mi. prior to reaching the outer marker. The 
same profile can also be synthesized by specifying a total of three-segments. 
By doing so, the altitudes at the ILS intercept and outer-marker points 
can be explicitly specified, thus ensuring at the outset that a level-flight 
segment of specified altitude and length will be included in the synthesized 
profile. The numerical values for four altitude waypoints are also shown in 
table 3. 


SYNTHESIZING SPEED PROFILE FOR CONTROLLED 

TIME OF ARRIVAL 

The speed profile for a specified controlled time of arrival is 
derived from interpolation between two other speed profiles, the maximum- 
speed profile and the minimum-speed profile. In the maximum-speed 
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profile, the aircraft is flying at maximum allowable speed at all times 
and the time-duration for the flight, t m in, is the mimimum achiev- 
able time of arrival. Similarly, in the minimum-speed profile, the 
aircraft is flying at minimum allowable speed at all times and the time- 
duration of the flight, t max , is the maximum time of arrival. Denote 
by td the desired time of arrival. If it is within t min and tmax, that 
is, 

tmin ^ td < tmax ( 15 ) 

then a speed profile can be synthesized to achieve the controlled time of 
arrival. 

From the perspective of air-traffic control it is desirable not only 
to maximize the range of arrival times but also to provide continuous 
control of them within the feasible range of times. This requires specifying 
a rule for generating a family of true-airspeed profiles which covers the 
entire speed envelope, from the minimum to the maximum speed of the 
aircraft. An appropriate parameterization of the speed profile can be 
developed by examining the airspeed envelope for a typical jet transport 
shown in figure 8. The figure indicates that the minimum true airspeed as 
a function of altitude lies along a constant-IAS contour. The maximum 
true airspeed as a function of altitude is determined both by maximum IAS 
and a maximum Mach number. The change from Mach number limited 
to IAS (or equivalently dynamic-pressure) limited flight occurs at about 
25,000 ft. The combined effect of these limits gives rise to the characteristic 
corner in the maximum airspeed envelope at about 25,000 ft. Typical airline 
procedures recognize this characteristic by recommending a constant-Mach- 
number descent at the higher altitudes (above about 25,000 ft) followed 
by a constant IAS descent at lower altitudes. Thus, an operationally 
and aerodynamically significant parameter for covering the speed envelope 
is a command value of IAS. The entire family of true-airspeed profiles, 
parameterized by the IAS is shown in the figure. If, at the initial altitude, 
the commanded IAS is larger than the maximum allowed at that altitude, 
the first part of the descent is flown at constant Mach number. The 
descent Mach number should be chosen near the operating limit in order 
to achieve a large variation in arrival time. For example, in figure 8 
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SPEED PROFILE 



Figure 8.— Constant CAS profiles. 


Mach 0.84 was selected. As the aircraft descends at this constant Mach 
number, IAS increases and will eventually reach the commanded value. 
The remainder of the descent is flown at the commanded IAS. The fastest 
possible descent first follows the Mach-number limit and then the IAS limit, 
thereby tracing out the maximum-speed envelope. At 10,000 ft altitude all 
profiles must decelerate to a maximum of 250 knots IAS (KIAS) to meet 
ATC rules. 

The parameterization of the speed profile as described above has an 
important effect on the method of computing the time-to-descent through 
a large altitude range. Since the true airspeed is not constant during des- 
cent, it is generally necessary to compute the time-to-descent by numerical 
integration method. Moreover, an iterative procedure, wherein each itera- 
tion step involves integrating a descent trajectory, must be used to com- 
pute the commanded indicated airspeed that results in a specified arrival 
time. 
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The complete speed profile for a realistic descent trajectory can 
consist of many segments of acceleration, deceleration, constant Mach 
number, and constant-IAS and constant- TAS flight. The method and al- 
gorithm described below provide for the synthesis of such complex speed 
profiles. 


Specification of a Complete Profile 

A speed profile is contructed from concatenation of speed seg- 
ments. A typical speed segment is shown in figure 9. It con- 
sists of three legs, typically an acceleration/deceleration leg, followed 
by a constant-IAS leg, and then a deceleration leg. The numeri- 
cal method used to synthesize a speed segment will be discussed 
later. 



Figure 9 — Three-leg speed segment-CIAS profile. 


A speed segment is specified by V' 0 the initial indicated airspeed, 
V' f the final indicated airspeed, V' c the constant indicated airspeed, and 
d 9 the length of the segment. 

Since synthesizing the speed segment involves integrating aircraft 
equations of motion, the indicated airspeed or Mach number must be 
converted into true airspeed before the integration. The corresponding true 
airspeeds are denoted V 0 , Vf, V cl , and V c2 . The latter two are, respectively, 
the true airspeed at the beginning and the end of the constant-indicated- 
airspeed leg. They may differ if the altitude changes while flying this 
leg. The true airspeeds are calculated from indicated airspeed by the 
relation, 
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" y l * 

1 

v = a(h)^ 

2 

7- 1 

( 

+ 1 ) 1 - 1 
V 2-7 P.1 T ) 

+ ‘) - 


(16) 


where a is the speed of sound, 7 is the ratio of specific heat, and P(h) is 
the atmospheric pressure (ref. 9). 

In terms of waypoints, a speed segment is specified by the coor- 
dinates of two adjacent speed waypoints, {sk,V' k ) and 
and constant indicated speed V’ c , where s* and s&. j-i are dis- 
tances from the touchdown point measured along the known horizon- 
tal path. The horizontal distance d 6 for a speed segment is there- 
fore 


ds — Sk s/t-j-i (IT) 

In general, the true airspeeds V c \ and V C 2 must be computed during 
integration since they are altitude-dependent and their corresponding al- 
titudes are not always known a priori. 

A specification for a typical speed profile consists of five speed 
waypoints, as described in table 4. 

Figure 10 shows the corresponding speed profile. The first seg- 
ment contains three legs and the remaining segments contain only 


TABLE 4.- SPEED WAYPOINTS FOR THE TERMINAL AREA 


Speed 

waypoint 

no. 

Distance 
from TD, 
n. mi. 

^ias 

Description 

Altitude 

Spl 

50 

300 

Cruise speed 

Above 10,000 ft 

Sp2 

40 

250 

Descent speed 

Above 1 0,000 ft 

Sp3 

30 

220 

Terminal-area speed 

Below 10,000 ft 

Sp4 

13 

170 

Approach speed 

Below 10,000 ft 

Sp5 

4 

120 

Landing speed 

Below 10,000 ft 
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Figure 10.— Typical speed profile. 


two legs. The length of the horizontal flightpath, which is also 
the distance from touchdown, is used as the independent variable. 
The horizontal distance locates the corresponding point on the syn- 
thesized horizontal flightpath and also that in the synthesized al- 
titude profile. The (x, y) coordinates in the horizontal flightpath and 
the altitude must be known at all times since the equation of mo- 
tion is altitude-dependent, bank-angle-dependent, and, possibly, wind- 
dependent. 

The coordinates of the first waypoint in the speed profile are 
(Spo , Vq ), the current aircraft distance from the touchdown point, and the 
current aircraft speed, respectively. The remaining waypoint coordinates 
are specified in table 4. 

Synthesizing a Speed Segment 

The three-leg speed segment illustrated in figure 9 is synthesized by 
(1) integrating forward from V 0 until V cl is reached; (2) integrating back- 
ward from Vf until V& is reached; (3) integrating forward from V c \ to V C 2 
to generate the constant IAS leg; and (4) computing the appropriate quan- 
tities such as the distance d Cl the time-duration, and the fuel consumption 
for the segment. 

Forward and backward integrations- Equations for integrating for- 
ward from V 0 at altitude h 0 to V cl at an initially unknown altitude 
h c i to generate the first leg in the segment will be described. The 
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TAS acceleration/deceleration is computed by using equations (18a) and 
(18b), 


V a = cos 7 — g sin 7 (18a) 

m 

for acceleration, and 

t'a = — — — — V w cos 7 — <7 sin 7 (186) 

for deceleration. In equations (18), T ma x and T m i n are the maximum and 
minimum thrust, respectively; D is the drag force; g is the acceleration of 
gravity; 7 is the flightpath angle; m is the mass of the aircraft, and V w is 
the component of the wind-rate vector along the heading direction. The 
wind vector V w can be written as the vector sum of its north (N) and east 
(E) components: 


V w = iV N + jV E (19) 

where V}y and V E may be functions of altitude. The time-derivative of V w 
is 



- dVj^r . - dV E l 

l — ~ — /l j— J —z~. — h 


dh 


dh 


( 20 ) 


If xp is the aircraft heading, a unit vector ip in the heading direction 
can be written as 


A. A A 

ip = i cos ip -\- j sin ip (21) 

o 

Then the component of V w along the heading direction is obtained by 
taking the dot product, 

V w =V W • ip = cos ip + sin ip)h (22) 

oh oh 
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Since the independent variable of the synthesized altitude profile and speed 
segment is the distance to go from touchdown measured along the syn- 
thesized horizontal path, it is necessary to rewrite equations (18a) and (18b) 
as a function of distance. Noting that 


y. 



(23) 


where V g is the ground speed; substituting equation (23) into equations (18a) 
and (18b) and integrating results in 


=/( : 


max 


— D V w cos 7 -+- g sin 7 


mV 0 


) 


ds -f Vo 


(24a) 


for acceleration, and 


V. = [( ~ ° - Y^212±l^l iL + (24t) 

j y mVg v g j 

for deceleration. Relationship between V a ,V g , and V w is expressed by the 
vector equation 


v f = v a + V w ( 25 ) 

The thrust is a function of engine pressure ratio (EPR), altitude, 
and Mach number. The maximum thrust is obtained by setting EPR 
equal to EPR ma x and, similarly, the minimum thrust is obtained by 
setting EPR equal to EPR m i n . The thrust model used in this simula- 
tion is identical to that used in reference 8, except for some minor 
simplifications. 
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The drag force D is computed by using equation (26), 


D=^p{h)SV a 2 C D (C L ,m,f r ,S. p ) (26) 

"where p{h) is the atmospheric air mass density, which is altitude-dependent; 
S is the wing reference area; and Cl is the drag coefficient. The drag 
coefficient is a function of the lift coefficient Cl, the Mach number, the 
flap deflection 6f, and the speed-brake angle. In the simulation, a drag 
model that approximates that of a Boeing 727 aircraft (as in ref. 8), was 
used. 

The lift coefficient Cl is computed by the relation, 

W 

Cl = — o < 27) 

^ p(h)SV g 2 cos 0 


where W is the weight of the aircraft and 0 is the bank angle. The 
value of bank angle is computed by the use of equation (1) for curved 
flight. 

The flap deflection 6j is 0° in accelerated flight, and is IAS- 
dependent in decelerated flight. The flap-deployment schedule, as a 
function of IAS, is approximated in figure 11. The flap deflection 
is 40* at or below 150 KIAS, and linearly decreases to 0° at 
210 KIAS. 
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Figure 11.— Flap deployment schedule. 
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If a segment requires decelerating and descending simultaneously, 
it is generally necessary to deploy speed brakes in order to achieve 
a reasonable deceleration rate. The effect of the speed brakes is to 
increase the value of the basic drag coefficient, Cd • The drag coefficient 
ACo t ssp owing to speed brakes is approximated as a constant value 
of 0.01 below the Mach number of 0.73. Its effectiveness decreases 
to zero at M = 0.95. Between these values, a polynomial in Mach 
number is used to calculate the drag increment. Speed brakes are 
deployed whenever the deceleration without speed brakes decreases below 
0.06 g’s. 

Equations (24a) and (24b) can now be rewritten for incremental 

steps, 


a,*+l 


f Tmax,i Di 

“ ^ mV 9}i 


Vw,% COSH + 9 sin 7* 




» + V^i (24 a’) 


for acceleration, and 


o.»+l 


f Di ^ w,i 

~ V 


cos it + 9 sin 


V a 


in n\ 

— U si + 


(246’) 


for deceleration. 

The ground speed in the denominator for the next integration step 
is updated by the relation, 


= \Va,i+l + 


(25’) 


The time-duration for traversing the horizontal-path length in each integra- 
tion step is expressed by the relation, 


A U 


= A Si/(^ 


2 cos 


) 


(28) 
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The incremental altitude is updated by the following relation, 



(29a) 


for level flight and by 


/i*+i = 


OWi -f V^sinq* 
2 


H- 


(296) 


for climb/descent. In equations (28) and (29), 7* is equal to 7 descent for 
descent and equal to iciimb for climb. 


O 

The fuel consumption / is the integrand / over time 



(30) 


Expressing equation (30) in terms of the horizontal distance as the indepen- 
dent variable yields 


/ = 



ds 


(31) 


o 

The fuel-consumption rate / can be expressed as a product of thrust T and 
thrust specific fuel consumption ( TSFC ), 


f — T (TSFC) (32) 

The value of TSFC depends on altitude, speed and temperature. It is based 
on the model developed in reference 8. The incremental fuel consumption 
A fi is derived by combining equations (31) and (32) and then expressing 
the result in incremental form, 
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Af i =(T{EPR,h i ,m i ), TSFC(h it m i: e i ))l~ ’! i+l f V,< As, (31’) 

where 0 is the atmospheric temperature and m is the Mach number. 
It will be used for computing the incremental fuel flow for both the 
acceleration leg and the deceleration leg. In the first case, thrust is 
computed by setting EPR equal to EPR max and in the second case, equal 
to EPRmin . 

Since the backward integration to generate the third leg in- 
volves deceleration only, a modified form of equation (18b) is used 
for the longitudinal equation. All other equations for backward in- 
tegration are essentially the same as those for forward integration 
with the exception that a minus sign is placed on the right-hand 
side of equation (18b) to reflect the negative time increment, that 
is, 

6 V a = — [ ^ mtw — — V w cos 7 — g sin 7 | |<5i| (186’) 

m 

Using distance as the independent variable yields 

— | cos 7 — g sin 7 | As 

AV a = — (186” ) 

Constant indicated airspeed leg- This leg is synthesized either as a 
gradually decelerating or constant- TAS leg depending on the altitude profile. 
If the altitude is decreasing, it is a decelerating- TAS leg, and if it is constant, 
it is a constant- TAS leg. The technique for synthesizing a decelerating- 
speed leg has already been discussed in forward integration. Specifically, 
for this case, use the deceleration equation and set the flap deflection to 
a constant value, since the indicated airspeed remains the same in the 
leg. 

By rearranging terms in equation (24), the explicit relation for the 
requisite thrust to maintain constant indicated airspeed is obtained, that 
is, 
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T { = m V gi — t+ * — - + V Wti cos ^ sin 7 , + D, 

As, 


(24c) 


The true airspeed V^+i, at the end of the integrating step is obtained 
by first determining the altitude h from the altitude profile from known 
integration step size As and then substituting the specified Vi as and 
appropriate cr(h) and P(h) into equation (16) to obtain the value of 


The technique to synthesize a constant TAS leg will be considered. 
The thrust for this leg is determined by balancing the longitudinal force 
equation 

T = D (33) 

since this case occurs only in level flight. Because the speed is constant, 
the fuel consumption, as expressed in equation (31), can be simplified and 
rewritten as 

( (TSFC)T is ^ {TSFC)Td c (84) 

J Vg Vg 

where d c is the horizontal distance subtended by the constant- speed leg. It 
is calculated from dj (see fig. 9), 


dc — ( ) d a df 


(35) 


In the presence of wind, the integral must be broken down into straight-flight 
and turn segments as dictated by the horizontal flightpath. Combining 
equations (34) and (35) and expressing the ground speed V g in terms 
of the airspeed and the wind speed by the use of equation (25) 
yields 


(TSFC) T (sk — s*+i — d 0 — df) 
V*-f 


(36) 


Constant Mach number leg- If the altitude profile is flown at con- 
stant altitude, then this leg degenerates to a constant- TAS leg. This is 
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usually not the case. The true airspeed of the aircraft V is a product of 
Mach number m and the speed of sound a, 

V = mo (37) 

Taking the time derivative of equation (37) and equating the time derivative 
of the Mach number to zero yields the relation 

t* (9cr f da 1[r . 

V a = = m ^r v sin 7 (38) 

u/i on 

For a constant 7, the thrust required to maintain a constant Mach number 
is expressed by the relation, 


Tmin = D + sin 7 (m — V — + W) (39) 

g on 

which is obtained by substituting equation (38) into equation (18) and 
solving for T m j n . At each integration step, the minimum value of thrust 
in equation (24b’) is replaced by the value calculated from equation (39). 
The quantity can be approximated by a constant value below the 

troposphere, typically 4.1 x 10~ 3 . Above the troposphere, fyfc — 0. 

Equations (28) and (29) are used for computing the incremental time 
and the altitude changes and equations (31) and (32) for computing fuel 
consumption. 

Summary of Algorithm for Synthesizing Speed Segment 

The algorithm can be divided into three steps, each of which is dis- 
cussed below. 

Step 1- From known V' k and fi 0) integrate forward until V' c is 
reached. First, compare V' k with V' c to ascertain whether it is a deceleration 
leg, a degenerate leg, or an acceleration leg, namely 

n > n (40a) 
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for a deceleration leg, 


( 406 ) 


for a degenerate leg, 


n = n 
n < n 


(40c) 


for an acceleration leg, where the initial speed V' k and the constant speed 
V' c are usually given as indicated airspeeds. The initial TAS V 0 is obtained 
from V' k by use of equation (16). Next set up the initial condition for 
forward integration by equating, 


v; = Vo 

(41a) 

h% — h 0 

(416) 

= iv^+vy 

(25”) 


From known distance s k determine the bank angle <p from the syn- 
thesized horizontal flightpath and the flightpath angle 7 from the syn- 
thesized altitude profile. Use equations (26) and (27) to compute the drag 
force D. 

If it is an acceleration leg, set EPR to EPR m ax to compute T ma x 
and use equation (24a’) to complete the integration step. If it is a decelera- 
tion leg, then set EPR equal to EPR m in to compute T m in and use equa- 
tion (24b’) to complete the integration step. After V di+i is known, up- 
date the altitude by using equation (29). Finally, compute the incremen- 
tal time and the incremental fuel consumption by using equations (28) 
and (31’), respectively. Continue the integration until the speed V' c is 
reached. 

If the altitude is still changing in the last step of integra- 
tion, then the last step size is determined by trial and error in or- 
der to complete the integration for the leg. Specifically, the TAS 
Vdi+i > s compared with an IAS V' e to terminate the integration. 
The constant speed V r c must first be converted into TAS before the 
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comparison, and the conversion is altitude-dependent. The altitude 
at the end of the leg, as expressed in equation (29), depends on 
Vfli+i. 

The time-duration for traversing this leg is obtained by summing the 
incremental time-durations, 


«1 = Af. (42) 

ft 

where i refers to the i th integration step. The distance subtended by this 
acceleration/deceleration leg is d 0 . It is the sum of all incremental integra- 
tion steps, 

d„ = A (43) 

ft 


Step 2- Similarly, from known V'j and hj , integrate backward until 
V' c is reached. The procedure for integrating backward is practically 
identical to that of integrating forward with three exceptions. First, 
it is not necessary to compare the final speed V'j fe+1 with the constant 
speed V' c> since the last leg is always a deceleration leg. Secondly, 
the incremental form of equation (24b’) is used for the integration, that 
is, 


Va i+l =- 


Tmin ~D V w cos 7 -J- g sin 7 


m V, 


\As -f Vai 


(24ft”) 


gi 


9 * 


where As, the integration step, is a positive quantity. Finally, a sign change 
in 7 is required in equation (29b) to reflect the backward integration, that 
is, 


^*+1 


(V gi+1 +V, i )sin(-7) , 
r r h, 


(29ft’) 


Denote by t$ and df the time-duration and the distance subtended for flying 
the third leg, then 

f 3 =^Ati (42’) 

i 
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= X>«i 

i 

Step 3- Compute the time-duration and fuel consumption for this 
segment. Denote by T* the time-duration for flying the k th segment in the 
speed profile, then 

Tk = t\ -j- t c t -h t c i -j- <3 (44) 

where t c t and t c i are the time-durations for flying the constant- TAS and 
constant-IAS legs, respectively. 

Denote by F* the fuel consumption for the k th speed segment and 
by fufc, and / 3 the fuel consumption in the first, constant-speed, and 
deceleration legs, respectively; then the total fuel consumption for the seg- 
ment is 

F k =/i +/*+/*+/* (45) 

The fuel consumption for the first or the third leg is obtained by sum- 
ming the incremental fuel consumption for each integration step. The 
incremental fuel consumption is determined by equation (3T). The 
fuel consumption for the constant-speed leg is computed by using 
equation (36). The distance subtended by this leg is determined by 
equation (35). 


Controlled Time of Arrival for Speed Segment 

The controlled time of arrival T oa is the sum of time-durations for 
traversing all speed segments in the speed profile, for example, 

T*a = £ h (46) 

k 

Each segment is specified by two waypoints, ( s*,V* ) and («*_(_! , ); 
some fixed quantities; and by the constant speed V' c , a quantity that 
can range between V' min and V' max . Denote by V' norn a selected 
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value of V' e ; then by judiciously choosing V' nom , a specified time for 
traversing the segment can be achieved. The value of V' nom is limited 
by 

n.in < V'„ m < V' max (47) 

Figure 12 shows the minimum, nominal, and maximum-time-of- 
arrival speed segments. The maximum speed segment consists of ac- 
celerating from V 0 to V max , flying at constant V max , and then decelerat- 
ing from Vmax to Vj) this is also known as the minimum-time seg- 
ment. The nominal-speed segment consists of acceleration/deceleration 
from V 0 to V norn , flying at a constant V nom , and then decelerating 
from V norn to Vf; it is also known as the nominal-time-of-arrival seg- 
ment. The minimum-speed segment consists of acceleration/deceleration 
from V 0 to V m i n , flying at constant V m * n , and then decelerating from 
V m i n to Vj ; it is also known as the maximum-time-of- arrival seg- 
ment. 



Figure 12.— Minimum, nominal, and maximum time of arrival speed segment. 


Denote by t nom , and t max the minimum, nominal, and 

maximum times of arrival, respectively. Also denote by td,k the 
desired time of arrival for the k th speed segment. Then synthesizing a 
specified time duration for the k th segment involves testing that 
is, 

fm»n — td,k — f max (48) 

and, if it is within bounds, proceed to select a speed V oat k such 
that 

||foo,fc (49) 
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where t oa ,A is the time-duration for flying along the segment specified by 
V oat k, and €t is the allowable time-error. 

An efficient technique for determining the speed V oai k to achieve 
the time-duration t oai will now be developed. From known ( V max , t m ), 
{Vnomytnom), and [V min , t max ), three points on the unknow n-speed- 
versus-time-duration curve can be located. These points and the 
curve representing the exact but unknown relationship are shown in 
figure 13. 



Figure 13.— Controlled time of arrival versus speed. 


An analytical approximation of the curve is expressed by the 
relation 


V* a> * — 


ci 


^oa,k 


+ 


C 2 

t2 oa,k 


+ 


C3 


^ oa,k 


(50) 


where t oa is the time duration for flying along the speed segment. 
The coefficients Ci,C 2 , and Cs are obtained by solving the matrix equa- 
tion 



( 51 ) 
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Taking the inverse of equation (51) yields 

c = W- 1 U~ l V (51’) 

where 



The ei's that appear in the elements of the matrices W and U are the 
reciprocals of time, that is, 


e, = (56a) 

* max 

e 2 = (566) 

* nom 

«3 = 7 J— (56a) 

*min 

The explicit solutions for the inverse of W and U are, respec- 


tively, 



(57) 
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M 


Tmax 

0 

0 


0 

P nom 

0 


) 

min ' 


■where A is the 3x3 co-minor matrix of W, 


( 6362(63 — e 2 ) eie 3 (ei— e 3 ) 

62 — 4 «| - e i 

63 — 62 61 — e 3 


62 61 (62 — 


c i 

62 


c 2 

6l 


) 


(58) 


(59) 


and d is the determinant of W ; 


i= e 3 «2(e3 — e 2 ) + eie 3 (e, — e 3 ) + e 2 ei(e 2 — ei) 


(60) 


After Ci, C 2 , and c 3 have been determined, an estimate V oa ,k of 
V oa> k can be computed by substituting td,k into equation (50). The actual 
value t oa ,k is obtained by synthesizing the speed profile using the estimated 
speed V oai k- Experience with this technique indicates that the accuracy of 
the arrival time determined in this first iteration estimate will be adequate. 
However, if the inequality in equation (49) is not satisfied, a second iteration 
based on computing a new set of coefficients, ci,C 2 , and c 3 must be 
carried out. This is done by discarding one of the four known pairs of 
values V oay k,t 0 a, k- We choose the three values of t 0Uy k that are clustered 
closest to tdk . Then the estimating and synthesizing steps are repeated as 
before. 


Synthesizing Controlled Time of Arrival Speed Profile 

The controlled-time-of-arrival speed profile is constructed from a con- 
catenation of synthesized controlled-time-of-arrival speed segments. The 
controlled time of arrival T oa for the entire speed profile is the sum of time 
durations for all segments, 


Toa=Y,‘ 


oa } k 


k 


(61) 
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The individual segment is specified by d k ,V maXik ,V nomjk , and V min<k , 
where d k is the horizontal distance subtended by the segment. The nominal 
speed for each segment is computed by the relation, 


Fnom.fc — Fmm.fc ) H - Vmin,k (62) 

The quantity X ranges from 0 to 1; a value of 0 means V nOTri)k takes 
the value of V m i nk and a value of 1 means it takes the value of 
Vmax,k • In order to minimize computation time, the desired time of 
arrival for all but the first segment is assigned a value of t n om,k, that 
is 

td,k = fnom.fc (63) 

corresponding to a value of X = 0.6. By this assumption, all 

but the first speed segment is generated in a single synthesizing 
step. The desired time of arrival for the first segment td, 1 is there- 
fore 

td, 1 — td — (td, 2 -f" ••• + td,k -f- — + td,nk ) (64) 

where nk is the total number of segments specified for the profile. If a larger 
range of arrival times than that provided by the first segment is required, 
then the second and subsequent segments will have to be included in the 
synthesis. 

As an example, the desired time of arrival for the second segment is 
computed by the relation, 


td, 2 = td — toa, 1 — (td, 3 ••• + td,nk) (65) 

Here, t oaj \ would be the upper or lower limit of arrival time. In 
addition to saving computation time, this restriction is also a judicious 
choice for two other reasons. First, the horizontal distance for the 
first segment is usually much greater than the sum of all horizontal 
distances for the remaining segments. In the example shown in table 1, 
the horizontal distance for the first segment is 103 n.mi., and for the 
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sum of the rest of the segments it is only 16 n.mi. Therefore, the first 
segment determines about 90% of the range of the controlled time of 
arrival. 


Second, at any time during the flight, a new desired time of arrival 
may he issued to accommodate the latest aircraft-arrival schedule which 
reflects the current air-trafflc situation. Consequently, it is necessary to 
resynthesize the overall speed profile at times. In this case, only part of the 
first segment has been flown and the subsequent speed segments remain the 
same. 


Figure 14 shows (solid lines) the three synthesized speed profiles 
involved in the controlled-time-of-arrival computation. They are the 
maximum-, nominal-, and minimum-speed profiles. Each profile is specified 
by a set of speed waypoints, in this case there are five, indicated in 
the figure by numbers and by a 3D trajectory. The speeds for the 


264 KIAS 



Figure 14.— Controlled time of arrival speed profiles. 


42 


maximum-speed profile are 264, 277, 260, 160 and 120 KIAS; for the 
minimum- speed profile they are 264, 226, 210, 160, and 120 KIAS. Their 
corresponding distances from touchdown are, respectively, 119, 28, 10, 
3, and 0 n.mi., respectively. The nominal-speed profile is specified by 
using a value of 0.6 for X. The synthesis is carried out so that speeds 
specified at a waypoint are not achieved until the aircraft arrives at the 
waypoint. 

Each profile consists of four segments. In the maximum-speed 
profile, the first segment consists of an acceleration leg, a constant- 
true- airspeed leg, a constant-indicated-airspeed leg, and a deceleration 
leg. Note that the constant-IAS legs appear as gradual decelerations in 
descent in figure 14. In the minimum-speed profile, the first segment 
consists of a deceleration leg, constant-true- airspeed leg, constant-indicated- 
airspeed leg, and a deceleration leg. All remaining segments consist of a 
constant-indicated-airspeed leg followed by a deceleration leg. The time- 
durations of flight for these two profiles are 1309 and 1529 sec, respectively, 
and for the nominal-speed profile the time-duration is 1384 sec. If the 
aircraft is to arrive in 1400 sec, the required initial indicated airspeed is 
253 KIAS (437 KTAS) or, equivalently, a value of 0.524 for X. These 
values were computed by using the technique described previously. A 
single application of equation (50) yielded an arrival-time error of less than 
0.5 sec. 


CLEARED FOR APPROACH COMMAND 

The purpose of this controller-initiated command is to simulate 
realistically the capturing of ILS localizer and glide slope by an unequipped 
aircraft. It is typically issued by the controller when an aircraft is on the 
base leg, between 8 and 30 n.mi. from touchdown. In response to this 
command, an algorithm must generate the horizontal, altitude, and speed 
profiles, including deceleration to landing speed at the appropriate point. 
Furthermore, the algorithm determines, before synthesizing a trajectory, 
if the initial position and heading meet the criteria for ELS localizer 
capture. 
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Synthesizing of the trajectory proceeds in the following four steps: 
(1) horizontal capture, (2) altitude capture, (3) speed capture, and (4) con- 
trolled time of arrival. In each case, there are criteria for testing the 
feasibility of capture. 


Horizontal Capture 

Figure 15(a) shows the geometry of the touchdown region in the 
horizontal plane. The points marked OM and TD are the outer marker 
and the touchdown points, respectively. The variables x, y are the 
coordinates of the aircraft position in the so-called outer marker coordinate 
system. It is defined by choosing the outer marker as the origin and 
the extended runway centerline as the x-axis. The angle ip is the 
aircraft heading, measured clockwise from the direction of the positive 
x-axis. 


X 



a) LOCALIZER CAPTURE REGION IN OUTER MARKER 
COORDINATES 



b) CONVENTIONAL CAPTURE HORIZONTAL FLIGHTPATH 
Figure 15.— Localizer capture. 
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Criteria for a feasible horizontal capture trajectory consist of (1) a 
preliminary position and heading test, (2) a conventional capture test and, if 
it is failed, (3) a TST test. The first test can be summarized in the following 
set of equations. 



x < 0 

(66) 


—20° < i>< 95° 

(67’) 

for y<0 and 

—95® < if)< 20° 

(67”) 

for y>0. 




Equation (66) requires that the aircraft location be in front 
of the outer marker when capture is initiated. Equations (67’) and 
(67”) define the allowable aircraft heading which depends on the sign 
of the y-coordinates, as shown in position o' and a”, respectively, 
in figure 15(a). The shaded area in the allowable heading region of 
figure 15 is a candidate for a conventional capture and the cross-hatched 
area is a candidate for a possible TST capture. This preliminary 
test, which requires no computation, gives necessary, but not sufficient, 
conditions for capture. After passing the first test, the next two 
criteria will determine the shape of the horizontal flightpath, if there is 
one. 


Figure 15(b) shows the shape of the horizontal flightpath for conven- 
tional ILS capture. It consists of flying straight with current aircraft head- 
ing to point P, then turning to merge without overshoot onto the extended 
runway centerline at Q. 

Criteria for a feasible conventional capture trajectory are sum- 
marized in equations (68) through (70): 

0 < %i <20 n.mi. (68) 
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0° < ip< 95°, y < 0 (69’) 

—95° < 1) < 0°, 2/ > 0 (69”) 

|V>I < *(*») (TO) 

In equation (68), z < is the x-axis intercept, which is defined as the point 
of intersection between the extended runway centerline and the x-intercept 
line, a directed line emanating from the aircraft position and pointing in 
the direction of aircraft heading. The criterion stipulates that this point 
must not have passed the outer marker and must be no farther than 
20 n.mi. from the outer marker. If it is farther than 20 n.mi., it is 
too early to exercise the cleared-for-approach command. Equations (69’) 
and (69”) define the shaded region in figure 15(a). Equation (70) limits 
the intercept angle as a function of the z-intercept coordinate in order 
to prevent unreasonable turns close to the outer marker. The value of 
intercept coordinate z* and the corresponding maximum allowable angle 
are tabulated in table 5. 

The tabulated value is based on the U. S. Standard Terminal 
Instrument Procedure with some modification (see ref. 10). Figure (16) 
shows the maximum allowed intercept angle as a function of the intercept 
distance to the outer marker, corresponding to table 5, for capturing the 
localizer. 


TABLE 5.- PERMISSIBLE CAPTURE REGIONS 


Intervals of intercept distance 
Xj from outer marker, n. mi. 

4' (max intercept angle), 
deg 

0 -0.8 

15 

.8- 1.0 

30 

1.0- 2.0 

45 

2.0 -3.0 

60 

3.0 -4.0 

75 

4.0 -5.0 

90 

5.0 and larger 

95 
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Figure 16.— Maximum intercept angle as a function of intercept distance to OM. 

The tabulated value limits the turning angle and thus the turning 
distance of the capture flightpath. For example, if the distance is 
between 0 and 0.8 n.mi., the angle of turn is limited to 15°. If the 
aircraft position and heading pass all the tests specified, a conventional 
capture trajectory is feasible, with one exception. The aircraft position 
on the line of intercept may lie between P and (see fig. 17). In 
this case, only a TST capture trajectory as illustrated in figure 1(c) is 
possible. 

Figure 17 shows the capture switching diagram for a typical ap- 
proach direction. If the distance between the current aircraft posi- 
tion and the intercept z* (the intersection of the current aircraft head- 
ing^ with the extended runway centerline) is greater than the distance 
Pxi (where P is the beginning of a turn for an STS pattern), a 
conventional capture is used. Beyond this point, only a TST cap- 
ture is feasible. At or beyond the intercept, no capture is pos- 
sible. 
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A singular case also occurs if the aircraft heading angle is exactly 

zero: 


if ip = 0; then x t = oo 


(71) 


In this case, the aircraft is flying parallel to the runway centerline and a 
TST type of capture must be used. The three waypoints specifying the 
trajectory are the aircraft position {x,y), the x-axis intercept (xj,0), and 
the outer marker (0, 0). 

If any of preceding tests is failed, then a TST capture may still 
be possible provided the aircraft position and heading pass a third set 
of criteria. Basically, these criteria require the initial aircraft position to 
be within a narrow cone centered on the x-axis, and the initial heading 
close to the runway direction. In practice, such initial conditions generally 
allow the aircraft to perform a successful ILS capture. But they can cause 
the test of the x-intercept x t to fall easily outside the range permitted by 
equation (68). 

Criteria for an allowed TST trajectory are summarized in equations (72) 
through (74): 
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1 < x < 30 n.mi. 


(72) 

(73) 

(74) 


\tp\ < 20° 

LOSmax < 20 ° 

Equation (72) stipulates that the aircraft ^-coordinate be at least 1 n.mi. 
from the outer marker and no more than 30 n.mi. from it. Equation (73) 
prevents the aircraft from making unreasonable turns and equation (74) 
confines the allowable aircraft position to within a maximum line-of- 
sight angle relative to the runway centerline. The line-of-sight angle 
is measured from the outer-marker location and is computed by the 
relation 

LOS = tan -1 ^ (75) 

x 

This trajectory is specified by the aircraft position and heading (x, y, ip), the 
coordinates of the outer marker, and the runway heading (0, 0, 0). At the 
outer marker the aircraft is lined up with the runway heading and proceeds 
to fly straight to the touchdown point. 


Vertical Capture 

Figure 18(a) shows a general altitude-capture profile consisting of 
three degenerate segments, each consisting of one leg, that is, (1) descending 
from current altitude to outer-marker altitude at a specified descent 
angle, (2) flying level to the outer marker for a specified distance 
di ; and (3) descending at a glide slope of 3° from outer marker to 
the touchdown point. The distance d is the distance subtended by the 
horizontal flightpath, the distance da is the distance subtended while 
descending at a specified flightpath angle 7 <*, and the distance d om> td 
is the horizontal distance between the outer marker and the touchdown 
point. 

The current aircraft altitude h a is compared with the altitude h v 
shown in figure 18(a) to determine the existence of an altitude-capture 
profile. It is computed by the relation, 
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( 76 ) 


hy — Ajfn "j - dd tan ^d 
dd == d di d orrii td 

where h orn is the altitude at the outer marker, and di is the minimum 
distance the aircraft must fly level in order to capture the glide slope. 
If the aircraft altitude is above h v then no capture is possible. If 
it is the same as h v then the altitude profile is the same as that 
shown in the figure. Finally, if it is below h v a level-flight segment of 
length 

Si, = h «- h ° (77) 

tan id 




b) SPECIAL CASE OF ALTITUDE PROFILE FOR 
GLIDESLOPE CAPTURE 

Figure 18.— Glideslope capture. 
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is added at the initial altitude h a before beginning the descent to altitude 
h om . A case of interest is di = 0, and id =3°. The altitude profile for 
this case is shown in figure 18(b). 

Speed Capture 

Figure 19 shows the speed (IAS) versus distance profile used for 
speed capture. The symbols S a / C , S a , and Si stand for current aircraft 
speed, approach speed, and landing speed, respectively. The distance d 8>rn i n 
is some specified minimum distance for decelerating from any allowable 
initial aircraft speed to approach speed before the aircraft has reached 
the outer marker. The criterion for the feasibility of a speed-capture 
profile is that d, the horizontal flighthpath distance, be greater than 
d s ,min + dom,TD- The length of the approach- speed leg depends on 
the current aircraft speed. The deceleration distances d a / c ,a and d a ,i are 
computed by forward integration after the speed-capture criterion has been 
satisfied. 
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Figure 19.— Nominal speed profile capture. 


Speed Profile for Controlled Time of Arrival 

The controlled-time-of-arrival algorithm previously described can 
be adapted for use in the near-terminal area in conjunction with the 
cleared-f or- approach command. This feature allows the controller to 
specify a landing time whenever he issues the cleared-for-approach com- 
mand. 
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The minimum', maximum-, and nominal-speed profiles for this 
case are shown in figure 20. The nominal profile, which is iden- 
tical to that in figure 19, is used when no time specification is 
given. If a time is specified, maximum-, minimum- and nominal- 

speed profiles are synthesized to determine the range of arrival times 
and to provide input data to the speed-iteration algorithm described 
previously. 



Figure 20.- Maximum and minimum speed profile for cleared for approach command. 

The algorithm for the synthesizing-speed profile, with or without 
specified controlled time of arrival, is shown in figure 21. If at any 
time a speed capture is not feasible, the air-traffic controller must vector 
the aircraft to a new position that allows speed-capture. The algorithm 
consists of the following steps: First, compare d with d 8l mi n + d 0 m,TD 
and, if it is greater than the sum, proceed; if not, then there exists no 
capture-speed profile. Second, synthesize the minimum- and maximum- 
speed profiles and compare the desired time of arrival td with f m »n 
and t m ax * (1) if td < tmi n , then the desired time of arrival is 
too short and again there exists no time capture; (2) if td > tm ax , 
the controller must increase the horizontal flightpath length by path- 
stretching to achieve the desired time of arrival; and (3) if neither of 
these limits is exceeded, then synthesize the controlled-time-of-arrival speed 
profile. 
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SYNTHESIZE 
NOMINAL SPEED 
PROFILE 
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SYNTHESIZE MIN., 
MAX. AND NOMINAL 
SPEED PROFILES 
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[(c) 

SYNTHESIZE 
CONTROLLED TIME OF 
ARRIVAL SPEED 
PROFILE 



< (b) 




Figure 21.— Logic for synthesizing speed profile. 
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Real Time Regeneration of Reference Trajectory 


In the real-time simulation the synthesized 4D trajectory deter- 
mines the three-dimensional position of the aircraft as a function of time. 
It locates the position of an aircraft target on the controller’s display. 
Since the simulation is conducted in real-time and the trajectory is syn- 
thesized in fast-time, a method for storing the trajectory and regenerat- 
ing it in real-time must be developed so as to produce realistic target mo- 
tion on the controller’s display. A brief description of the method fol- 
lows. 

The horizontal and vertical trajectory is stored in a table listing 
the coordinates of the beginnings and ends of turns and their turning 
radii, as well as the length of each turn and straight-line segment. 
The vertical trajectory is specified as a function of distance along the 
horizontal path as previously explained. The synthesized speed profile 
determines the motion of the simulated aircraft target along the 3D path. 
During fast-time synthesis, the airspeed profile is stored in an array of 
coefficients as follows. In the TAS leg, the speed is constant and is 
invariant with time. In the constant-IAS legs, the speed is expressed 
as a linear function of time, and in the acceleration, deceleration, or 
constant-Mach-number legs, it is expressed as a quadratic of time. The 
coefficients of the quadratic are obtained by a least- square- fit technique, 
with the initial and final points of the leg weighted heavily to preserve the 
continuity between legs in the speed profile. Equations for these legs are, 
respectively, 


V = constant 


(78) 


V — €L\ * 4 " #2 t 


(T9) 


V — bi -f- b 2 f + b 3 t 2 


(80) 


where "constant,” ai , a 2 , b x , b 2 , and b 3 are the coefficients. 
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The synthesized speed profile is also used to determine the 
time for the beginning and ending of turns in the horizontal trajec- 
tory and the beginning and ending of altitude changes in the altitude 
profile. 

ENERGY-RATE MODEL FOR AIRCRAFT SIMULATION 

The ability of the aircraft to change speed and altitude simul- 
taneously is limited by the minimum- and maximum-energy rates, 
which are functions of thrust and drag. Although the synthesized 
trajectory takes into consideration the limits of energy rate, con- 
trolling a simulated aircraft to fly along the synthesized trajectory 
usually requires an energy rate that is different from that in the 
reference trajectory. For example, the winds, atmospheric condi- 
tion, and the performance model may be different from that assumed 
in the synthesis. If the requisite energy rate is outside the limits, 
then some adjustment of rate of speed change or flightpath angle or 
both is necessary. This section describes a method for handling this 
problem. 

The total energy of the aircraft is the sum of the kinetic energy and 
potential energy, 


The energy rate is obtained by taking the derivative of equation (81), 


The altitude rate can also be expressed in terms of V a and flightpath angle 


E = -?- + h 


(81) 



(82) 


g 


h—V a siB7 


(83) 
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Combining equations (82) and (83) yields 


* V a v a . 

E = (- V sin 7 

9 


( 84 ) 


Normalize equation (83) by dividing it by V a forms the dimensionless 
quantity 

k V a . . /otx 

— +sm7 (85) 

Va 9 

The left side of equation (85) "will be referred to as the normalized energy 
rate, that is, 

k n = — + sin 7 (86) 

9 

Taking the wind into consideration, equation (86) becomes 


b K , . , dV w cos 7 

E n b sin 7 + — - 

9 at g 


(87) 


Combining equations (19) and (87) and rearranging terms yields, 


V a . 

b sin 7 = 

9 


L max 


-D 


W 


(88a) 


Y + si ^= Tm 'w ° ( 88i > 

The left sides of equations (88a) and (88b) are, respectively, the maximum 
and minimum normalized rate, that is, 

£»,»», = (89 a) 
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K,min = - 7 ’ - m H 7 - (m 

These two energy rates form the upper and lower bounds of k n , 

&n,min ^ k n ^ k n,max ( 90 ) 

If b n is within bounds, then the requisite speed and altitude changes 
called for in their respective profiles can be met and the simulated 
aircraft can follow the synthesized trajectory exactly. Otherwise, 
some allocation of energy rate for speed and altitude change must be 
made. 


Since it is more important to minimize altitude error than speed 
error (to clear obstacles or ensure altitude separation), tracking the 
altitude profile is given priority over tracking the speed profile. Thus, if 


either 

E n > E njmax 

(91a) 

or 

E n < E njTn i n 

(916) 


then some recomputation of V a is necessary. If the requisite energy rate 
is greater than the maximum energy rate the aircraft can deliver, then 
set 

E n == E n ,max (92a) 

and adjust V a to reflect this change: 

V a = g (E n>max — sin 7 ) — V w cos 7 (926) 

Similarly, if the requisite energy rate is less than £ ni mtn, then 
set 

kn ^ k n ,min (93a) 
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I 


( 936 ) 


and again adjust V a to reflect this change: 

== 9 {&n,min sin'y) 

Equations (87) and (88) are used to test the energy rate when controlling 
an aircraft to fly along a synthesized 4D trajectory. If it is outside the 
limits, then equation (93) is used to adjust the rate-of-speed change V a . The 
readjustment of the speed profile is limited by the maximum and minimum 
speeds of the aircraft. The thrust is adjusted to maintain the aircraft 
within its speed envelope. The result will be a time-of-arrival error which 
may necessitate intervention by the controller if it creates a scheduling 
conflict. 


CONCLUDING REMARKS 


A set of algorithms for synthesizing 4D fuel-efficient trajectories has 
been developed for simulating an advanced air-traffic-control environment 
in a terminal area. The algorithms include realistic modeling of aircraft 
performance and fuel consumption. 

Flexibility in synthesizing trajectories is achieved by defining three 
types of waypoints. Waypoints in the plane specify the horizontal profile; 
altitude and speed waypoints are defined with respect to distance from 
touchdown, measured along the horizontal profile. This approach allows 
the horizontal profile of an approach route to be modified without affecting 
the altitude and speed profiles. 

The major emphasis in this report is on the development of the 
time-control algorithm. Time-control is obtained primarily through the ap- 
propriate choice of the speed profile. First, the maximum range in touch- 
down time is computed by integrating the point-mass equation of motion 
along the upper and lower speed boundaries. Second, if the desired ar- 
rival time falls within this range, the appropriate speed profile is computed 
by an iterative procedure. The computationally intensive method of in- 
tegrating equations of motion to generate speed profiles is necessary in or- 
der to model realistically the effect of aircraft performance and altitude- 
dependent winds. The algorithm applies to a wide range of operational 
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conditions. Thus, specification of arrival time can be given as early as a few 
hundred nautical miles from touchdown or as late as a few tens of nautical 
miles. 


Although the algorithms described herein have been developed 
specifically for simulation of advanced air-traffic-control systems, they can 
also be applied to the solution of a variety of on-board trajectory-synthesis 
problems. 
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